Phonons and d-wave pairing in the two-dimensional Hubbard model 
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We analyze the influence of phonons on the d I .2_ !/ 2-pairing instability in the Hubbard model on the 
two-dimensional square lattice at weak to moderate interaction U, using a functional renormalization 
group scheme with frequency-dependent interaction vertices. As measured by the pairing scale, the 
B lg buckling mode enhances the pairing, while other phonon modes decrease the pairing. When 
various phonon modes are included together, the net effect on the scale is small. However, in 
situations where d-wave superconductivity and other tendencies, e.g. antiferromagnetism, are closely 
competing, the combined effect of different phonons may be able to tip the balance towards pairing. 
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The two-dimensional Hubbard model is one of the 
most-studied models in context with high-temperature 
superconductivity in the layered cuprates. An important 
question is whether electronic interactions of the Hub- 
bard model alone can provide a sufficient pairing strength 
explaining the high critical temperatures and large en- 
ergy gaps observed experimentally. Furthermore, in par- 
ticular for large values of the Hubbard interactions, there 
is a strong competition between various ordered states 
such that one may want to have an argument why super- 
conductivity prevails in a large parameter region for most 
of the materials. Although the lattice degrees of freedom 
were thought to be irrelevant for the high-T c problem for 
a long timeQ, and phononic signatures in the electronic 
properties are still debated intensively 2], an additional 
phononic contribution to the pairing seems to be a nat- 
ural way to enhance the superconducting pairing against 
other competing electronic correlations. 

A theoretical analysis of this question in the Hubbard 
model at large values of the onsite interaction is diffi- 
cult. The impact of phonons on the pairing interaction 
has been addressed in various ways with partially con- 
tradicting results[3, Here we propose to analyze 
the situation at weak to moderate values of the coupling 
constant. This will allow us to obtain some qualitative in- 
sights. The method we use is the functional renormaliza- 
tion group (fRG) . Previously, this method has been used 
to classify the leading instabilities of the weakly coupled 
Hubbard model without phonons 
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band fillings when the Fermi surface (FS) is not nested, 
a d x 2_ y 2-w&ve superconducting instability is obtained in 
a large parameter window. The main driving force for 
these superconducting tendencies are antiferromagnetic 
(AF) spin fluctuations. Here we add phonon-mediated 
interactions to the bare Hamiltonian. We analyze the 
changes in the critical energy scale for the Cooper insta- 
bility and in the competition with other states. 

The fRG scheme we use is an approximation to an ex- 
act flow equation for the one-particle irreducible vertex 



functions of a many-fermion system^]]]- The quadratic 
part of the fermionic action is supplemented with a cut- 
off function which restricts the functional integral over 
the fermions to the modes with dispersion |£(fc)| > A. 
For the 2D square lattice we use a ^'-parameterization, 



-2t(cosk x + cosky) — At' cos k x cos k y — \x with 



nearest and next-nearest neighbor hoppings t and t' and 
chemical potential \i. The fRG flow is generated by low- 
ering the RG scale A from an initial value Ao ~ band- 
width. Thereby momentum shells with energy distance 
A to the FS are integrated out. In the approximation 
we use, the change of the interaction vertex is given by 
one-loop particle-hole (including vertex corrections and 
screening) and particle-particle pairs where one interme- 
diate particle is at the RG scale A while the second one 
has j£(fc)| > A. Higher loop contributions are generated 
by the integration of the RG flow. The flow of the self 
energy will be analyzed in a later publication (see also 
an d as in previous studies its feedback on the flow 
of the interaction vertex is neglected. 
The interaction vertex Va(&1j &2, ^3) depends on the gen- 
eralized wavevectors of two incoming particles (k\ and 
£2) and one outgoing (£3) particle with wavevector, Mat- 
subara frequencies and spin projection ki = (fcj,u;,, Sj). 
The dependence of this function on three wave- vectors is 
discretized in the so-called A^-patch scheme, introduced 
in this context by Zanchi and Schulz 7]. This scheme 
takes advantage of the tested fact that for standard Fermi 
liquid instabilities the leading flow is rendered correctly 
by projecting the wavevectors k\,ki and k^ on the FS 
and keeping the variation of V\{kf , fcf, fcf, • ■ • ) when 
the k[ are varied around the FS. Hence, one calculates 
V\(ki, ki, £3, . . . ) for k±,k2,k3 on the FS and treats it 
as piecewise constant when k\,ki, and kj, move within 
elongated patches stretching from the origin of the Bril- 
louin zone (BZ) to the (±7t, ±7r)-points. In order to treat 
retarded interactions we have to go beyond the previ- 
ous works which neglected the frequency dependence |l3|]. 
We divide the Matsubara frequency axis into M sections. 
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The aim is to approximate the decay of a phonon propa- 
gator above a characteristic frequency ujq. Below we show 
results for 32 BZ patches and M = 10. The minimal fre- 
quency spacing is ujq = 0.2t. The frequencies for which 
the vertices are computed range between ±5ujq. The fre- 
quencies of the dispersion-less phonons considered here 
are taken to be less than ujq. We have checked that other 
reasonable choices do not change our qualitative findings. 

The RG flow is started at an initial scale A with ini- 
tial interaction V\(ki, k 2 , &3j W\, oj 2 , W3). What is typi- 
cally encountered at low T is a flow to strong coupling, 
where for a certain flow parameter A c one or several com- 
ponents of Va(A:i, k 2 , ^3, wi, u>2, UJ3) become large. At this 
point the approximations break down, and the flow has to 
be stopped. Physical information about the low-energy 
state is obtained by analyzing which coupling functions 
and susceptibilities grow most strongly. For standard 
Cooper instabilities, the critical scale A c at T = is pro- 
portional to the critical temperature T c . 

For pure Hubbard interactions, the initial vertex at 
scale A is VF (hi, k 2 , ^3) = U. For phonon-mediated 
interactions, we add a retarded part, leading to 



under 90 degree rotations, 



V Ao (k 1 ,k 2 ,k 3 ) = u-J2 



gi{ki,k3)g l {k 2 ,k i )uj lfi 
(wi - w 3 ) 2 + uh 



(1) 



Here, <?i(fci, k$) is the electron-phonon interaction for an 
(Einstein) phonon mode i with frequency ujq scattering 
an electron from k\ to k s . Obviously, if the product of 
the g's in the numerator does not produce sign changes, 
the main effect of the phonon part is to reduce the ef- 
fective onsite interaction U . The coupling strength of 
the mode can be measured by the FS average Xi — 
2 Eg,* ) \9i& fcO| 2 /KiVol Efe Moti- 

vated by current issues in the high-T c problem, we an- 
alyze various different phonon modes, idealized as dis- 
persionless. First we consider a Holstein phonon with a 
fc-independent coupling ga ]stein(k,k') = g. As a next 
step we analyze a coupling which only depends on the 
transferred wave- vector q — k — k' , 



|3Buck(g1| 2 = <?Lck(cOS 2 ^ 
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A coupling of this type was used by Bulut and 
Scalapino|l4| in their analysis of the out-of-plane mo- 
tion of the planar oxygens (buckling mode) in the lan- 
guage of a one-band Hubbard model. They and vari- 
ous other authors 0, 0, 0] pointed out that the sup- 
pression of this coupling for large momentum transfers 
q ~ (it, 7r) leads to an attractive d-wave Cooper pair- 
ing potential. The same mode was also considered using 
multiband models involving the planar or full copper- 
oxide structure [13, Then the out-of-phase c-axis mo- 
tion of the planar oxygens on x- and y-bonds, known as 
Big-mode, gives rise to a coupling with a sign change 



SBi g (fc,fc') = SBig v(l + cosfc^Xl + cosk' x 



— + cos k y )(l + cos k'y) 



(3) 



roughly 60meV, and the coupling stren gth is discussed 
to be as high as A apcx = 0.5 and higher [ 



In optimally doped cuprates this mode shows up at 
~ 36meV and its coupling strength is taken as Aeig = 
0.23[I3. In-plane breathing modes have also been ob- 
served and discussed in the cuprates. Here the pla- 
nar oxygen atoms next to a copper site move to- 
wards or away from the copper atom. One obtains 
a coupling (again with q = k - k') ^breathe® | 2 = 
Sbroatho ( sin2 + sin2 if) E3- Tne frequency for this 
mode in optimally doped samples is taken as 70meV, 
and the FS average of the coupling strength is quoted 
as Abrcathc = 0.02^^. Finally there is a c-axis vi- 
bration of the apex oxygen above or below the copper 
with a coupling g apcx (fc, k') = g apcx (cosfc K -cosk y ) ■ 
(cos k' x — cos k' y ) . The bare g ap ex(k, k') does not sup- 
port an attractive d-wave pairing component, but it has 
been argued^] that screening can generate a d-wave at- 
traction. The frequency of this mode in the cuprates is 

i gth is 

Drl9l. Below we 
use the conversion 50mcV= Q.lt. 

First we describe the flow without phonons. We choose 
a curved FS near the van Hove points with n = 0.83 par- 
ticles per site. Note that the doping dependence at weak 
coupling is quite different from the behavior at large U, 
so no strong conclusions about the doping dependence in 
the cuprates can be drawn. There are two main effects of 
the interaction: tendencies towards AF spin-density wave 
(SDW) and towards d-wave pairing. The FS with the 32 
discretization points, and the flows of d-wave pairing and 
AF-SDW susceptibility Xdw and \af are shown in Fig. 
n As in previous studies without frequency dependence, 
for these parameters Xdw grows most strongly toward low 
scales, but one clearly observes tight competition with 
antiferromagnetism. Therefore, for these situation with 
the FS near the saddle points, alternative interpretations 
of this multi-channel instability have been considered0. 
Reducing U or the nesting makes the c?-wave pairing more 
dominant. In the right panel of Fig. Q] we display the fre- 
quency dependence of \dw which shows the build-up of 
a zero-frequency peak when the instability is approached 
at low scales. We also plot the FS-averaged d-wave pair- 
ing interaction with zero total incoming frequency versus 
the frequency transfer. The peak is quite broad, signaling 
a large pairing 'Debye' energy scale ujd ~ t. It roughly 
tracks the frequency dependence of \af at wavevector 
(tv,tt). A similar behavior has been found in dynamical 
cluster approximation (DC A) calculations 22J. 

Now we include phonon modes. The critical scales A c 
for the (i-wave pairing instability are shown in Fig. [21 
a) as function of the dimensionless coupling strengths A 
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FIG. 1: (color online), a) RG flow of the susceptibilities 
Xdw(w) (thick (dashed) line for u) = (cj = 0.5i)) and 
XAF(<f, <j-" = 0) at wavevector q = (7r,7r) (thin solid line) for 
t' = -0.25t, U = 2.5t and /i = -O.Mt with AT = 32 FS 
points, T — O.Oli. The marks on the A-axis denote where the 
largest coupling reaches 5t and Wt. b) Frequency dependence 
of Xdsc(<^) at the scale where Kiax reaches 32t (thick solid 
Une), Wt and 5i (dashed lines). The thin solid line shows the 
rescaled xaf(w) at V ma x = 32t. The open circles show he FS- 
averaged d-wave pair scattering V£ w (u), —U),w m ) versus trans- 
fered frequency u> (with outgoing frequency u m = ±0.1t). 

for the various modes included separately at U = it. 
The phonon-mediated interaction has two effects which 
compete. First, a momentum-dependent structure can 
develop which can generate a d-wave component in the 
pair scattering, and enhance (i-wave pairing. Second, the 
attractive part of the retarded interaction reduces the 
effective onsite repulsion |2l|. which can disfavor spin- 
fluctuation-induced pairing. The competition between 
these can be illustrated by the Holstein coupling. 

At least for U < 6t it is known that the Hol- 
stein coupling g (k, k') is suppressed by the electronic 
correlations |3l t|- The suppression is strongest at large 
k — k' w (7r,7r)- The fRG for a Holstein mode with fre- 
quency 50meV reproduces this trend 20]. In principle, 
this generates a d-wave component in the pair scattering 
and A c for the d-wave pairing instability should increase 
by adding the Holstein phonon. However, the fRG finds 
a reduction of A c . The reason is the suppression of the 
initial Hubbard interaction by the Holstein phonon which 
outweighs the additional d-wave attraction. 

Next we consider the buckling phonon @ with 
frequency 50meV. Now, already the initial phonon- 
mediated interaction is attractive for d-wave pairing, as 
the scattering with k — k'w (tt, tt) is more repulsive than 
for k — k' pa 0. But again, the reduction of the effective 
initial repulsion is too strong, and the net A c is lower 
than for pure electronic interactions. A c vs. A behaves 
similarly to the Holstein case and is not shown in Fig. [21 
We conclude that it is not justified to simply add phonon- 
mediated pairing interactions on top of unchanged spin- 
fluctuation-induced interactions in a BCS gap equation. 
The interaction between these two channels needs to be 
considered. Similar trends are found for breathing and 
the apex oxygen modes. For stronger coupling to the 
breathing mode, Abrcathc > 0.5, the d-wave pairing in- 



FIG. 2: (color online), a) Critical scales A c for the d- 
wave pairing instability vs. electron-phonon coupling A, for 
T = O.Oli, fj, = -0.95t, t' = -t/4 and U = 3t. Holstein 
mode: diamonds), breathing mode: squares, Si 9 -mode (Eq. 
[3J: crosses, and apical mode: circles, b) Critical scales A c 
for the d-wave instability vs. U. The dashed line is without 
phonons. Crosses (squares) for only the B\ g (apical) mode 
included, triangles for the breathing, B\ g and apical mode 
together. All data for T = O.Oli, fi = -0.95t, t' = -t/4. 

stability gives way to a s-charge-density wave instability 
with modulation wavevector (tt, tt). The only phonon we 
studied with a positive effect on A c for d-wave pairing 
is the Big buckling mode. Since gBig(k, k') vanishes for 
k — k' w (tt, 7t), it supports d-wave pairing while it does 
not suppress the onsite-£7 due to its d-wave form factor. 
In fact, the fRG energy scale for an AF instability at 
perfect nesting and half filling is enhanced by the B\ g - 
phonon. In addition, the bare gmg gets enhanced by 
the Hubbard interactions 0, |2^| during the flow. Hence 
the single £?i ff -mode added to the Hubbard model can 
increase the energy scale for d-wave pairing considerably. 
For n = 0.83, U = 2.5t and A Blg = 0.23[l8|, the increase 
is as high as 55%, for U = 3t only 16%. With breath- 
ing, Big and apical mode included together (with the As 
and frequencies cited above, Fig|2]b)), A c is still slightly 
increased compared to the case without phonons. How- 
ever, such a small increase may be affected by details or 
uncertainty about the As for different phonons. 

Next we analyze the competition of the d-wave super- 
conductor with antiferromagnetism. As shown in Fig. ^ 
a), the pairing susceptibility Xdw overtakes the growth 
of the AF spin susceptibility xaf only very close to the 
instability. When the -B lg -phonon is included 10 a)), A c 
grows. Also xaf is increased at a given scale A, but 
Xdw dominates more clearly. This means that the Bi g 
phonon is indeed beneficial for d-wave pairing correla- 
tions. In addition, the coupling to the i?i g -mode also in- 
creases the d-wave charge fluctuation tendencies at small 
wavevectorsjl^. This could lead to an additional break- 
ing of the fourfold symmetry of the FS0,|2j|. This would 
not suppress the pairing instability altogether, but could 
reduce A c by pushing density of states away from the FS. 
If we now add the breathing and the apical mode (Fig. 
El c )), Xaf and A c get reduced again. However, Xdw is 
less affected by the weak breathing mode and the apical 
mode which barely changes the d-wave pairing. Hence for 
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FIG. 3: (color online). RG flow of Xdw (solid line), xaf 
(dashed) and d-wave charge susceptibility at small q (dashed- 
dotted) for U = 2.5t, t' = -0.25t, fi = -0.941 a) with only 
the Big-mode, Asi 9 = 0.23, b) with also apical and breathing 
modes included, A S i s = 0.23, \ ape x = 0.5, A brcathc = 0.02. 



U = 2.5t, compared to the case without phonons or with 
the Big phonon alone, Xdw dominates even more clearly. 
Similar trends are seen at U = 3t. Hence, at least at 
weak coupling, there is the possibility to enhance <i-wave 
pairing in energy scale and with respect to competing in- 
stabilities by coupling to the right mix of phonons. 
Finally, commenting on possible caveats for our analysis, 
we note that selfenergy corrections have been neglected. 
The reduction of the quasiparticle weight at small ener- 
gies could potentially reduce the pairing strength even 
further 15]. DCA results^ for U = 8t and Holstein 
phonons are consistent with this. The fRG results give an 
upper bound for the energy scale for d-wave pairing. Note 
that at large U near half filling the coupling to Holstein 
phonons may actually increase the AF susceptibility |2^| 
contrary to our weak coupling results. The reason is 
the formation of a heavy polaronic quasiparticle which 
is more effective in an already correlation-narrowed band 
at large U. Nevertheless the decrease of the pairing scale 
seems a common feature at weak and strong coupling. 

In conclusion, we have analyzed the influence of var- 
ious phonon modes on the d-wave pairing instability in 
the 2D Hubbard model at weak to moderate coupling. 
Most phonons studied reduce the energy scale for the 
instability of the Fermi liquid state by reducing the ef- 
fective onsite repulsion. This effect outweighs possible 
enhancements of the d-wave pairing scale due to the 
wavevector dependence of the electron-phonon coupling. 
Spin-fluctuation- and phonon-mediated pairing interac- 
tions are not additive. The only mode studied here which 
enhances the energy scale for the d-wave instability is the 
Big buckling mode. Due to its d-wave-type wavevector 
dependence it does not suppress the local onsite repulsion 
and therefore does not harm the spin-fluctuation mech- 
anism. For moderate U ~ 2.5t and average coupling 
Asi g = 0.23, the increase of the pairing energy scale is 
more than 50%. This increase is reduced when other 
phonon modes are included. Notably, for the parameters 
used here, where without phonons the d-wave pairing and 
AFM were in very tight competition, the combined effect 



of three phonons is able to establish the dominance of d- 
wave pairing relative to other instabilities. 
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